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I. 


Acconplishineint of Technical Objectives 


A. Develop digital filter programs using large pixel windows 

This capability has been con^letely developed and implemented. 
The BAFl and BAF2 programs are quite versatile and effective. Complete 
specifications are contained in the software description. 

B, Develop geometric correction programs necessary for registration 
of X-ray images 

This objective has not been met to date. It is anticipated 
that this program will be iir^ileniented during the continuing effort on 
the N-Ray /X-Ray Study. 

C. Develop an operational image processing system capable of handling 
large image arrays 

This capability has been achieved using an approach which 
analyzes sub-arrays of the entire array. Upon receipt of equipment 
from the N-ray/X-ray study, this capability will be extended farther 
and optimized on even larger arrays; without additional funding. 

D, Develop more efficient data transfer routines to be used by 
existing image processing systems 

This objective has been attained and realized by every program 
listed in the software description. Execution times of one-tenth (0.1) 
the previous execution time are typical. 
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II. Overall Assessment 

Capability for X-ray enhancement at the KSC Data Analysis 
Facility has been significantly increased as a result of this study. 
Information, heretofore mattainable, can now be extracted from 
X^ray imagery used in KSC support, A repertoire of software to 
optimally analyze various X-ray imagery has been successfully 
developed. Further X-ray enhancement capability (e.g. geometric 
correction, remoting of processed image to user site) will exist 
after completion of the N-ray/X-ray study. X-ray imagery used in 
KSC launch and support activities can new be analyzed more rapidly, 
thoroughly and correctly. 
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III. X-r-ay Enhancement l^cedures Manual 


Goal I : 

Tool: 
Method : 


Example: 


Sharpen edges to inprove visibility of inage detail. Extract 
subtle edges which are obscured by gradual shading effects. 
Minimize the effects of some shading nonlinearities. 

High-pass Filter 

1. Digitize X-ray using Video Averaging to increase signal- 
to-noise ratio in the image. 

2. Using the BAF2 digital filter, input an appropriate 
high-pass filter transfer function — suppressing the low 
frequencies . 

3. Using Contrast Stretch (Piecewise Linear or Linear) re- 
distribute the gray levels in the filtered image for 
optimum presentation to the eye. 

4. If excessive Vgraininess" exists in output image, use 

the Variable. Low-pass Filt^ program to slightly attenuate 
the highest spatial frequencies. 

In an analysis of an X-ray of the Space Shuttle solid rocket 
motor, the High-pass Filter and High-Eirphasis Filter success- 
fully revealed more definitive information about defects. 
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Goal II; Deblunr an X-ray to recover edge information and detail. 

Tool: High-Bnphasis Filter 

Method; 1. Digitize X-my using Video averaging to increase signal- 
to-noise ratio in image. 

2, Using the BAF2 digital filter, input an appropriate 
hi^-emphasis filter transfer function — boosting the 
hi^ frequencies on the original image. Extent of de- 
blurring is determined by: Ca) the slope of the transfer 

function, first, and then (b) the size of the filter'. 

Gain of the transfer function should be just below the 
saturation level. 

Exanple: In an X-ray analysis of a pmrp, the Hi^-Enphasis Filter 

successfully deblurred the X-ray. Ifere detail was visible 
after the restoration. The initial blurr function was caused 
by an out-of -focus condition which resulted when the pump 
was not in the focal plane of the X-ray system. 
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Goal III : Remove random noise (e.g. graininess or salt and pepper 
effect) frcm an image to reveal information previously 


masked by noise. 

Tool : Low-pass Filter 

Method: 1. Digitize X-ray using Video averaging to minimize any 

random noise from the camera/digitizing system. 

2. Using the BAF2 digital filter or Variable Low-pass 
Filter, input an appropriate low-pass filter transfer 
function — suppressing the hi^ frequencies. Cutoff of 
the transfer function should be no lower than necessary 
to just reveal features of interest — further filtering 
will begin blurring the image unnecessarily. 

Example: When enlarging a small area of many X-rays, film grain in- 

troduced random noise to the image. The Low-pass Filter 
successfully minimized this effect and revealed information 
previously undetectable. One particular exaanple involved 
small cracks in a crane, hook X-ray . 
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Goal IV; Extract only the edge information from an image to see if 
bounderies of features are correct. Display the different 
degrees of edges. 

Tool: Nine Point Gradient 

Method: 1. Digitize X-ray \jising Video averaging to increase signal- 

to-noise ratio in the image. 

2. Using Linear Constrast Stretch, distribute image gray 

, V * 

levels to fill the dynamic' range of solid state memory 
(256 gray levels) . 

3. Apply the Nine Point Gradient program. Output gray levels 
are proportional to gradient (edge) magnitude of input 
imge. 

Exanple: In an analysis of an X-ray of a printed circuit board, the 

Nine Point Gradient successfifLly detected defects in the 
board and electrical components. 
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Goal V: 


To display small gray level vaipiations using false colcsrs 
for increased differentiation. To extract areas of similar 
gray levels using colors to represent each pattern. 

Tool : Density Slicer 

Method: 1. Digitize X-ray using Video Averaging to increase signal- 

to-noise ratio in the image. 

. 2. Using Density Slicer in the automatic mode, set a cursor 
over the area of interest and specify the number of lini- 
formly spaced intervals to be displayed in false colors. 
In the nanual mode, input the gray level bounds fcsr each 
interval. Each output false color represents one finite 
range of gray levels. 

Example: In an analysis of metal test bar X-rays for a corrosion study, 
the Density Slicer was used to map the extent of corrosion 
growth quickly and effectively. 
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Goal VI: To improve the contrast of an X-ray, or any particular feature 


within an X-ray — yielding more information about the content 
of the X-ray. That is, to improve the display of small changes 
in gray levels so that they may be better differentiated from 
their surroundings. 

Tool: Contrast Stretch (Linear and Piecewise Linear) 

Method : 1 . Digitize X-ray using Video' Averaging to increase signal- 

to-noise ratio in the image. 

2. In Linear mode, place cursor over area of interest, or 
input gray level bounds. The input gray level range will 
be expanded to fill the full dynamic range of memory 
(256 gray levels). In Piebewise Linear mode, enter break- 
point coordinates of the input/output gray level transfer 
function. Output gray levels are distributed according 
to the piecewise linear transfer function. 

Example: As previously mentioned, the contrast stretch programs were 
used as a post-processing enhancement for most all X-ray 
enhancement investigations'. Nearly every X-ray can benefit 
from improved contrast. 
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Goal VII: Minimize the effect of shading nonunifomity within an 
X-ray 

Tool: Shading correction, Ranp Mask,. Background Mask, Linear 

Combiner 

Method: 1. Digitize X-ray using Video Averaging to increase 

signal-to-noise ratio in the inage. 

2. If a background mask is available, apply Shading 
Correction after the nask is registered, 

3 i If a mask is not available, use Ramp Mask (for a linear 

distortion) or Background -Mask (for a nonlinear distortion) ' 
to generate a mask containing the nonuniformity in shading. 
For additive distortions, use Linear Conibiner to subtract 
the mask from the original — effectively correcting the 
nonuniformity. For multiplicative distortions, use 
Shading Correction to correct the original by computing 
the ratio of the original to the mask. 

Exanple; In a corrosion analysis of a Space Shuttle Thermal Protection 
System tile X-ray, the Ramp ffesk and Linear Combiner programs 
effectively corrected a severe 'shading uniformity caused by the 
X-ray scanner system. 
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Goal VIII; Output enhanced X-rays on a high-re'solution hard copy 


Tool: 
Method : 


E>aiiple: 


(film or half-tone image on line printer paper) . 

Gould Print (halftone) , Micro-D Tape Write (film) 

1. Use the Gould Print program to dump the enhanced image 
on the line printer. Single or double page resolution 
may be selected. Contrast Stretch may be necessary for 
best output on paper. 

2. If film output is desired, use the Micro-D Tape Write 
program to dump the enhanced image onto magnetic tape. 
The tape may be input to the film printer of the micro- 
densitometer. A very high resolution transparency will 
be generated. 

These hard copy programs were used successfully on many 

X-rays . 
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Goal IX: 

Tool: 

Method: 


Example: 


Input precisely digitized X-rays from a microdensitometer 
for a more detailed analysis . 

Micro-D Tape Read 

1. Use microdensitometer to generate a magnetic tape' con- 
taining the digitized X-ray. Quantization is very detailed: 
256 gray level resolution and 12.5 micron spot size capa- 
bility. 

2. Use Micro-D Tape Read program to read the magnetic tape 
information into solid-state memory for further processing. 
X-Y geometric corrections may be performed. 

The Micro-D Tape Read program can be used for any X-ray of 
228 X 228mm size. Improved spatial frequency quantization 
and gray level quantization will be realized. 
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IV. Brief Software Description 

A. Previously existing software used in X-ray erihancement support 

1. Densitometer 

Purpose - To display the "density" (gray level) distribution 
of a chosen line for a given channel, field number, and 
raster line. 

2. Video Scanner Averaging 

Pilose - To reduce the noise (improve signal-to-noise ratio) 
caused ty random processes in the television camera scanning 
system. 

3 . Digital Theme Print 

Purpose - To reproduce a selected binary theme on the line 
print er/plotter . 

4. Micro-D Tape Read/Write 

Purpose - To input digitized image data from -the hi^-reso- 
lution drum-scanning microdensitometer — enabling further 
processing by computer. Dumps data in Image 100 memory in 
a format readable by the microdensitometer — enabling film 
printing of the memory image. 

5. Tape Registration 

Purpose - Register a video image to another image or reference — 
enabling further processing of images on same spatial reference. 

6. Shading Correction 

Purpose - Correction of illumination nonuniformity over the 
field of view for better radiometric accuracy. 
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7. Single Cell Signature Acquisition 

Purpose - Determine the single cell signature from the statis- 
tioa.l computations (frequency distribution-histograms) of the 
spectral intensity distribution within the training site. 

J 

8. Single Cell Histogram Display 

Purpose - Display histogram of single cell signature acquisition 
(# picture elements vs. gray level) allows modification of gray 
level limits. 

9. Window Program 

Purpose - Expand (enlarge) a selected area by a factor of 
•2 to 9 and 'train on the area. 

10. Ratio Mode ! 

Purpose - Reduce multiplicative errors in images. 

f 

11. Contrast Stretch 

Purpose - jfedmi.ze data range over ayailable memory bits — for 
visual enhancement of data for viewing. Small gray level 
variations are made large. 

12. Density Slicer 

Purpose - Examine the spectral limits for a given image and 
slice the range of gray level limits • into a specified number 
of slices. 
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B-. Software written as a result of the X-ray Enhancement Software 
Developnent and Test (attached) 

1. Nine Point Gradient 

Purpose - To display the approximate nagnitude of the gradient 
(rate of change in gray level) of an image for edge distri- 
bution analysis. 

2. Laplacian 

Purpose - To display the Laplacian operator (second order 
derivative) of an inBge for edge enhancement — approximates 
hi^-pass filter over small areas. 

3. Variable Low-pass Filter 

Purpose - To provide a variable degree of smoothing over an 
image by suppressing high frequency noise. 

4. Piecewise Linear Transfer Characteristic (Contrast Stretch) 
Purpose - To display the contrast of gray levels in an image 
in the optimum nanner, by allowing the user to construct any 
desired transfer characteristic. 

5. Gould Print of Full Screen 

Purpose - To print a halftone reproduction of any image stored 
in memory. 

■ 6. Background Mask Generator 

Purpose - To aid in correction of illimiination nonuniformities 
by generating a mask representing the distortion. 
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7 . Ramp Mask Generator 

I 

Purpose - To aid in correction of illumination nonuniformities 
by generating a ramp mask representing the piecewise linear 
distortion. 

8 . Linear Combiner 

Purpose - To allow synthesis of any linear combination of 
images for image enhancement. 

9 . BAF Spatial Filter 

Purpose - To allow flexible and efficient digital filtering 
for enhancement of an image's spatial frequency components. 
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Nine Point Gradient 


PROGRAM SGQRGES: GRAD. FTN, GRAD. MAC 


ORIGINAL PAGE IS 
OP POOR QUALITY 


PURl’OSii: "niis program approximates the magnitude of the gradient of an image 
stored in one channel of video memory and outputs it to the same or another 
channel. 


INPUTS: The user inputs the channel of the image for which the gradient is com- 
puted as well as the channel w'hich receives the output image. 

OUTPUTS: The computed gradient is output to the selected channel of video memory. 

PR OCESSING: Each element Vi, k of the output Jmage matrix V is computed 
according to the following algorithm. (Matrix f represents the image matrix for 
which the gradient magnitude is approximated. ) 

The first quantity approximates the magnitude of the vertical partial derivative 
and the second quantity that of the horizontal partial derivative. To reduce compu - 
tation time, these two magnitudes are summed rather than computing the 
Pythagorean distance of the two. 

It should -be observed that the V matrix is shifted diagonally one pixel with respect to 
the input image matrix. 
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Nine Point Laplacian 

.PROGRAM SOURCES: LAPL. FTN, LAPL. MAC 

PURPOSE: This program approximates the Laplacian operation on an image stored 
in one channel of video memory. The computed Laplacian image is placed in a 
selected channel of video memory. 

INPUTS: The user enters the channel of the image for which the Laplacian is 
computed as well as the channel which receives the output image. Also, he 
enters an integer gain factor between 0 and 7 inclusive which scales the output 
by that power of two. 

OUTPUTS: The computed Laplacian is output to the selected channel of video 
memory. 

-*2 

PROCESSING: Letjjs represent the output image by a matrix V and the input 
image by a matrix f. The matrix is computed according to the following 
algorithm : 

— <<>■». ” 3 

N is a user inputted scale factor between 0 and 7, inclusive. A bias of 128 is used 
so that positive and negative values of the quantity in brackets can be accommodated 
by the video memory. 

It should be observed that the V matrix is shifted diagonally one pixel with respect 
to the input image matrix. 
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Variable lowpass Filter 


PFOGEAM SOUECES: VM?LOW.FIN, VABLCW.MAC 

PUKPOSE: This program provides a variable degree of smoothing of an image. 

INPUTS; The user ii^uts the charnel of the image to be smoothed as well as the 
channel nuriber on which to output the snoothed image. Mso, an integer between 
0 and 100 inclusive is input which r^resents the percentage of smoothing desired. 

COTPOTS: The smoothed data is output to the selected channel of video matiory. 

PECXZEISSING; Let X3S represent the_output image by a matrix S and the input 
image by a matrix f . The matrix S is ocnputed according to the following 
algorithm: 


N is the user irputted percentage of smoothing desired. * 

Observe that the S matrix is .shifted diagonally one pixel with respect to the 
irput image matrix ?T 








160 



18 



origin ALPA^® 

OF POOR 


PROCJRAM NAMI-:: Pn<:CKWlSii UNIiAR I'RANSI-IiR CflARACTERISTIC 

PUKPOSI'): This |:)r<)griim nllows the user lo construct any desired transfer 
charactcrisl 1C and apply it to any image in the video memory. 


PROGRAM COMPONENTS: PLTC. FTN, PLTC. MAC 

INPUTS: The user enters ordered pairs of numbers (integers between 0 and 255, 
inclusive) which correspond to points in the cartesian plane describing input 
vs. output gray levels. (Entry of a number out of the 0 to 255 range is interpreted 
as signalling the end of input. ) Also, the user specifies the input channel to which' 
the transfer characteristic is applied and the output channel of memory in which 
the result is placed. 

OUTPUTS: The transfer characteristic as specified is drawn on the Tektronix 
CRT terminal. The image resulting from application of the transfer characteristic 
to the input image is placed in the selected channel of video memory. 


PROCESSING: To describe the mechanism of the transfer characteristic, consider 
the following plot on a cartesian plane. The X axis represents the input gray level 
and the Y axis the output gray level. Moving in the direction of ascending values . 
of X, the user specifi^ points are connected with line segments to form the 
transfer characteristic. The output image is formed by replacing the pixel value 
specified on the X axis with the value specified on the Y axis. 
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PROGRAM NAME: GOULD PRINT OF FULL SCREEN 

PURPOSE: This program prints a halftone reproduction of any of the video channels 
on the GOULD 5000 graphics line printer. Two sizes of output are available. 
Halftones are represented in a 4 by 4 dor matrix. 

QOMPONENTS: GOUl’RT. FTN, GOUl'R’l'. MAC, GOULD. MAC 
(GOULD. MAC is a G. E. Image 100 supplied driver for the GOULD-5000. ) 

INPUTS - The user inputs the desired size of the output, single sheet or double sheet 
Also, he enters the channel of video to be printed. 

OUTPUTS - The specified channel of video is printed on the GOULD 5000 line printer 
in either one or two sheet format, depending on the selected output size. 

PROCESSING: Two stages of processing occur: formatting the video data to match 
the specified output size; conversion of the 8 bit gray level to the corresponding 
4 by 4 dot matrix used to represent a halftone. 

In the case of the single sheet output, the '512 by 512 pixel image matrix is reduced 
to 256 by 256 pixels by averaging clusters of four pixels. For the double sheet 
output, the 512 by 512 pixel matrix is split into two 256 by 512 pixel matrices, ' 
corresponding to the left and right halves of the image. Each of these halves is 
subsequently processed and output. 

Each element of the 256 pixel line is converted to. a 4 by 4 dot matrix. Sixteen 
distinct dot patterns are used to produce sixteen distinct gray levels on the output. 
The upper four bits of the eight bit pixel value select which of the sixteen dot 
patterns is to be output. Each 256 pixel line thus produces four 1024 element lines 
on the GOULD printer formatted such that each 4 by 4 cluster of dots represents 
a halftone. 
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PROGRAM NAME: BACKGROUND MASK GENERATOR 

PURPOSE: Generates a mask by extending a horizontal or vertical profile of an 
image across an entire channel of video display memory. 

COMPONEN'PS: BGDMSK.Ii'PN 

INPUTS: The user selects the input channel' from which the mask is to be gen- 
erated and the .output channel on which the mask is to be placed. Also, he selects 
either a vertical or horizontal profile and positions the cursor to designate where 
the profile is to be taken. 

OUTPUTS-: The profile is extended and a mask generated on the specified output 
channel. 

PROCESSING: A horizontal (vertical) line is read, the position of which is specified 
by the cursor. On the specified output channel this line is printed in every horizontal 
(vertical) line on the screen so that the mask is generated. 
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1'KUC;kaM NAMI':; KAMI’ mask GI<:NERA'rOK 

PURPOSE: Extracts a ramp function from a selected portion of a horizontal 
or vertical line of data and generates a mask by extending this ramp profile 
across the entire screen. 

COMPONENTS: RMPMSK. FTN 

INPUTS: The user selects the input and output channels and either a horizontal 
or vertical mask. In addition, he positions the cursor (crosshair mode) over the 
portion of the line (vertical or horizontal) that he wishes to have the ramp function 
generated from. 

OUTPUTS: The calculated mask is output to the selected channel. 


PROCESSING: 


The ramp function is calculated as follows: 

y; s (XU^\) -for 

Y*- = 




where y* is the i*"^ pixel of the computed ramp function. (Values for ^‘greater 
than 255 are set equal to 255. ) ^The slope of the ramp, a, is computed using 
least squares linear regression according to the following algorithm: 



^ W* ( * 

e>~ (xt 


where: m and n are the boundary pixel numbers of the 
cursored portion of the selected input line. 


y* is the i^^ pixel of the input line. 




This algorithm for the slope a "is derived from the general formula for fitting a 
straight line to a set of points using the least squares method; 


Ol 


^ a. 
2L /)Ca 
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Once computed, this ramp function is printed in every horizontal or vertical line 
on the output channel (depending on whether a horizontal or vertical mask has been 
selected) so that the ramp mask is generated. 
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PROGRAM NAME: ONE TO FOUR CHANNEL LINEAR COMBINER 


DESCRIPTION: This program allows the user to construct an image which is a linea 
combination of up to four other images. The user specifies how the combination is 
to be done. 


PROGRAM COMPONENTS: LINCOM. FTN, LINCOM.MAC 
INPUTS: The program requests the following information 

(1) The number of channels to combine (1 to 4). 

(2) Which channels to combine. (Any combination of the five video 
channels.) 

(3) The weight for each channel. (A real number between -1. 0 and 1. 0, 
inclusive. ) 

- (4) Gray level bias. (An integer between 0 and 255, inclusive. ) 

(5) Output channel. (Any of the :(ive video channels. ) 

OUTPUTS; The combination specified is computed and output to the selected channel 
of the video memory. 

PROCESSING; The output image is computed as follows; 




A 

where: fo is the output image matrix 

N is the ^i^mber of input channels combined 

Ec is the 1 ” input image matrix ^ 

g; is the assigned weight corresponding to the image fi 

b is a 512 by 512 matrix,- each element of which is equal to 
the specified gray level bias. 
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1 . 


Introduction 


The BAF spatial filter program is supplied in two versions, 
BAFl and BAF2, The two versions differ only in the mechanism- for inputting- 
filter characteristics, as described in. section 1. 2. 2. 2. 

1. 1 Filter Type 


The BAF spatial filter operates as a convolution type filter. The 
operation is that of a running weighted summation. 

Fig. 1 demonstrates the operation of the running weighted sum- 
mation. Assume we have selected the summation to be done over a 3 by 3 
pixel window. The first part of Fig. 1 represents this window by a rectangle 
enclosing a 3 by 3 pixel area of the original image. The second part of the 
figure represents the weighted summation, where each pixel in the window 
is multiplied by the corresponding element of the filter weight matrix and 
these products summed, producing the output pixel. The last part of Fig. 1 
shows the position of this pixel in the output image. Each pixel of the output 
image is computed in this way, that is, centering the window over the corres- 
ponding pixel of the input image and performing the weighted summation. 

1. 2 Filter Design 


The filter weight matrix is designed by the user for each run 
of the filter. The user specifies' the size of the matrix (and the correspond- 
ing window) and supplies the specifications for the filter weights. 

1.2.1 Filter Size 

The filter matrix is a square matrix of dimension N. The value 
of N is selected by the user and can be any odd integer between 3 and 19, 
inclusive. Thus, the maximum filter size is a 19 by 19 filter weight matrix. 

1.2.2 Filter Weights 

Radial symmetry about the center filter weight is imposed so as 
to prevent phase shift in the output image. This symmetry is illustrated in 
Fig. 2 for a 5 by 5 matrix, in which weights with the- same number are forced 
to be equal. The user has the choice of inputting the spatial filter weights 
directly or inputting weights describing the 'frequency response of the filter. 

1.2. 2.1 Spatial Filter Weights ■ 

In version BAFl, the user can assign values- to the filter weights 
directly. Fig. 3 illustrates the assignment procedure for a 5 by 5 filter matrix. 


26 



The program will request the weights by the X and Y coordinates, which 
correspond to the numbered pixels in the figure. (The values entered must 
include a decimal point. ) The remainder of the filter weight matrix will be 
assigned values according to the symmetry described above. 

1.2. 2. 2 Frequency Domain Input 


Rather than assign the filter weights directly, the user can 
describe the frequency response desired and the appropriate filter matrix 
will be computed. Each of the two versions, BAFl and BAF2, provides its 
own method of specifying the frequency response. 

1.2. 2. 2.1 BAFl Frequency Domain Description 


In BAFl, the user specifies the elements of the G matrix defined 
^ the transform pair (3. 4) of the appendix Imposing radial symmetry in the 
G matrix ensures radial symmetry in the gjnatrix, as is desired. Fig.' 4 . 
illustrates the assignment of values to the G matrix. The program will 
request the values by the horizontal and vertical frequency components, which 
correspond to the elements numbered in the figure. (Decimal point must be 
included in input values. ) The remainder of the Gjnatrix will be assigned 
values according to the symmetry of Fig. 2. The g matrix will be computed 
according to the transform (3. 4) of the appendix, resulting in a filter matrix 
with the frequency response as specified. 

1.2. 2. 2.2 BAF2 Frequency Response Description 


Often one is concerned only with radial frequency response. 

BAF2 enables the user to specify a radial frequency res^nse curve. The 
program will sample this curve to obtain values for the G matrix which will 
then be processed as described for BAFl above. , , , 

The curve produced for the radial frequency response is seg- 
mented. The user specifies the endpoints of the segments. Fig. 5 gives an 
example. In this case, the user has specified three points. To form the 
response curve, consecutive points are joined with line segments. (The 
program orders the points according to the x-axis value, i. radial frequency. ) 
The program will sample this curve to obtain values for the G matrix. * 


* Given u, the horizonta l frequency, and v, the vertical frequencv, the radial 
frequency is 
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Fr equ ency Domain Characteristics 




The relationship between a given filter matrix g and its spatial 
frequency response G is developed in the appendix and restated here: 




. V i 1 


for -M ^ a, b, k, n ^ M and 


N = 2M+1 


Each element ^ corresponds to the frequency response of the filter at the 
frequency ), where N is the filter size (matrix dimensions) and T is the 

sampling periM of the digitized image. We have imposed via radial symmetry 
that 


Observe that ^/N corresponds to the horizontal digital frequency in cycles per 
sample and to the vertical digital frequency. 

2. 1 Filter Size Considerations 


The choice of size of the square filter matrix to be used has con- 
siderable effect on the type of frequency response curve that can be obtained. 

The number of elements in the frequency response raa^ix G is 
equal to the number of elements in the corresponding fulter matrix g. Thus for 
a filter dimension of N there are elements in each matrix. However, radial 
symmetry allows us to specify only 1/2 (M+l)(M+2) unique values in either matrix, 
where M = . Increasing N thus increases the number of points we can specify 

in the frequency response curve of the filter. 

The interval between adjacent points in the G frequency response 
matrix is ^/NT. In terms of digital frequency this is 1/N cycles per sample. 

The highest digital frequency that can be specified horizontally or vertically is 
M/N 

W • 

1 9 

For N = 3 this frequency is -^/3; for N = 19, this frequency is -fg- - 0.474. As N 
is made larger, this upper frequency specification approaches 1/2. This is in 
keeping with the Nyquist criteria which places the-highest digital frequency com- 
ponent of a properly sampled signal at 1/2. 
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3. 


Method of Solution 


Fig. 6. illustrates the procedure used in deriving the filter and . 
performing the operation. 

The user enters the size of the filter desired as well as the input 
and output channels of video memory. If the filter is to be described in the 
frequency domain, the G matrix is input, _either directly in BAFl or via the 
curve sampling technique of BAF2. The G matrix is transformedjo the g 
matrix. If the filter is to be described in the spatial domain, the g matrix is 
input directly. 

At this point, the floating point matrix g is scaled by an appropriate 
power of two to bring ilk into a range that can be represented with maximum pre- 
cision by a 16 bit integjer array. The scale factor is saved for later restoration 
and the scaled matrix g is fix^ to 16 bits. A bias value is computed which when 
summed to the output of the filter operation places the expected value of the output 
at 128 so as to be compatible with the 8 bit video display memory. The filter 
operation is implemented with the fixed 16 bit representation of g. 32 bit fixed 
arithmetic is used for the intermediate multiplications and summations so that 
no precision is lost in this stage. (Each pixel value is represented by 8 bits 
and the corresponding filter weight by 16 bits. ) The fixed point Output values 
from this operation are shifted right or left the proper number of times to com- 
pensate the scale factor. These values are each summed to the bias value to 
produce the values that are displayed. Negative display values are forced to 0; 
values greater than 255 are forced to ^5. All output values are thus compatible . 
with the video display memory. 
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APPENDIX 

Mathenatical developenent of BAF spatial filter 
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